Correlated theoretical uncertainties for the one-jet inclusive cross section 
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We discuss the correlated systematic theoretical uncertainties that may be ascribed to the next- 
to-leading order QCD theory used to predict the one-jet inclusive cross section in hadron collisions. 
We estimate the magnitude of these errors as functions of the jet transverse momentum and rapidity. 
The total theoretical error is decomposed into a set of functions of transverse momentum and rapidity 
that give a model for statistically independent contributions to the error. This representation can 
be used to include the systematic theoretical errors in fits to the experimental data. 
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Predictions of the Standard Model are typically made 
with the aid of next-to-leading order (NLO) perturbative 
calculations (or sometimes with NNLO calculations) . Ev- 
idently, these predictions are not exactly equal to what 
one should measure if the Standard Model is correct. If 
we have an NLO calculation, we leave out NNLO and 
N 3 L0 contributions, etc. We also leave out contributions 
that are suppressed by a power of the large momentum 
scale of the problem. Of course, we do not know exactly 
how big these contributions are: if we could calculate 
them, we would include them in the prediction. Never- 
theless, we can estimate the size of the corrections. They 
then constitute "theory errors" in the prediction, which 
are quite similar to experimental systematic errors in the 
measurement. 

In this paper we distinguish between errors associated 
with higher order contributions and power suppressed 
contributions to the cross section, which we call theory 
errors, and errors associated with our imperfect knowl- 
edge of the parton distribution functions needed for the 
prediction. Estimated theory errors are needed in two 
contexts. First, if an experiment does not agree with the 
theoretical prediction within the experimental statistical 
and systematic errors, then we need to see if there is 
agreement within the combined experimental and theory 
errors and the errors from the parton distributions used 
in the prediction. In the case that the disagreement is 
outside of the combined errors, then we have a signal for 
new physics. 

The second context in which we need estimated the- 
ory errors is in the determination of parton distribution 
functions from experimental measurements. The theory 
errors give a contribution to the errors that we associate 
with the parton distribution functions that emerge from 
a fit to the data. Evidently, if we do not include the- 
ory errors, the resulting errors in the parton distribution 
functions will be too small. Additionally, if for one kind 
of process the theory errors are large while for another 
kind of process the theory errors are small, then we will 
give the large-error process too much weight in the fit. 

In this paper, we provide an estimate of the theory 
error for the one jet inclusive cross section d 2 a/dPTdy 
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in hadron-hadron collisions, where Pt is the transverse 
momentum or "transverse energy" of the jet. and y is the 
rapidity of the jet. There is good data for this process 
from the CDF and DO experiments at Fermilab, including 
careful estimates of the experimental systematic errors. 
Estimates of the theory errors are needed to accompany 
the estimates of the experimental systematic errors. 

We warn that there is no unique method to estimate 
theory errors. Thus our task is to provide a method that 
is defensible if not necessarily optimal. We seek to pro- 
vide an estimate in a form that includes the correlations 
from one {Pt,u} point to another. 



II. GENERAL SETUP 

We treat theory errors in a fashion that is similar to 
that used for correlated systematic errors in the exper- 
imental results. We use next-to-leading order quantum 
chromodynamics (QCD) theory to make predictions for 
the one-jet inclusive cross section 1 
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In the calculation, one uses Monte Carlo integration so 
that there is a random statistical error for each point 
{Px,y}- We do not include these statistical errors in 
the analysis here since they are typically quite small (say 
2%) and one can reduce them by running the program 
for a longer time. If we wished to include the errors from 
fluctuations in the Monte Carlo integrations, that task 
would be straightforward because the statistical nature 
of these fluctuations is known. 

We will start our investigation by studying jet produc- 
tion corresponding to the Tevatron Run 2, with y/s = 
1960 GeV, as a function of Pt and y. We will display the 
results for y = {0,l,2}as functions of Pt; we also present 
formulas for the Pt and y dependence, from which esti- 
mated errors for the specific kinematic ranges used by 
CDF and DO can be inferred. 

We need estimated errors that can be used in a statis- 
tical analysis. However, we do not have at hand a statis- 
tical ensemble of worlds in which terms beyond those in- 
cluded in the NLO theory vary. Thus we make estimates 
that we hope are reasonable but that can and should be 
subject to debate. 

We formulate the treatment of theory errors as follows. 



1 Specifically, we use the program of Ref. 1 , although there are 
other programs that can give the same results. The code is avail- 
able at |http:/ /zebu. uoregori.edu/~sop er/EKSJets/jet. html 



We let 
da 



dP T dy 



da 



dP T dy 



NLO 



l + Y^XjfAP^y)] 



(i) 



Here the functions fj(Pr,y) are definite functions, while 
the Aj are unknown parameters. Thus Aj/j(Pr,y) rep- 
resents an unknown theoretical contribution that might 
modify the NLO theory. We treat the Aj as Gaussian 
random variables with variance 1. That is, the size of 
the uncertainty with label J is represented by how big 
fj{Pr,y) is. If one thinks of this as representing an 
imaginary ensemble of worlds in which theory calcula- 
tions come out differently, then these worlds all have the 
same fj but the Aj vary. 

We will propose to use just a few functions fj. We 
offer the following defense of this strategy. Consider a 
simplified case of a cross section that is a function of 
just one variable, Pt- If we were to believe that the 
uncertainty in the prediction of this cross section is of 
order, say, 10%, but we have no idea of what the shape 
of the true cross section is within a 10% band about the 
prediction, then we would choose many functions fj(Pr), 
each of size 0.10, but with each being non-zero only in 
a very tiny range of Pt- This approach is illustrated in 
Figure [l]-a) ; such a view seems to us unreasonable. 

Experience with various perturbative and non- 
perturbative contributions teaches that they are smooth 
functions of the relevant variables, Pt in this case. This 
arguably more reasonable scenario is illustrated in Fig- 
ure[l]-b). As illustrated by the three curves, 2 one contri- 
bution beyond NLO could be flat, amounting to a con- 
stant "K factor," another might be a smoothly increasing 
function of Pt, while yet another might be positive at 
high and low Pt and negative in between. However, we 
judge it unlikely that a currently uncalculated contribu- 
tion contribution would have multiple maxima between 
low and high Pt- 

Thus we seek a few functions fj(Pr, y) that have some 
dependence on {Pt, y} and represent, as best we can de- 
termine, our understanding of the character of uncalcu- 
lated contributions. In the following sections, we analyze 
several sources of theory errors and associate them with 
functions fj(P T ,y). 



III. PERTURBATIVE UNCERTAINTY 

The main source of uncertainty at large jet transverse 
momentum, at least in our estimation, is the fact that 



2 Specifically, in this figure we use the functions 
/i(P T ) = 0.1, fr{P T ) = 0.08 logCFr/Af), and f 3 (P T ) = 
0.06 |[log(P T /Af)] 2 - 0.1 j where M= 150 GeV. These curves 
are for illustrative purposes only, and the fj(Pr) functions 
differ from the set fj(PT,y) we will use to parameterize the 
correlated systematic uncertainties. 
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P T (GeV) P T (GeV) 

Figure 1: illustration of (a) uncorrelated and (b) correlated theoretical errors. In (a), the total error is about 10% for all Pt, 
but the error at any Pt is not correlated with the error at nearby points. In (b), there are just three functions fj(Pr) giving, 
again, about a 10% total error at any one Pt- Because the fj(Pr) are smooth functions, the theoretical error at a given Pt 
will be smoothly related to the error at other Pt values. 



we have calculated only at NLO, leaving contributions 
from higher orders of perturbation uncalculated. We es- 
timate this uncertainty using the dependence of the com- 
puted cross section on the renormalization and factoriza- 
tion scales. We present this estimate in this section. In 
the following section, we check this estimate using an 
independent method involving threshold effects. 



A. Error estimate from scale dependence 

The first ingredient in our estimation of theory errors 
is based on the traditional method in which one evalu- 
ates the dependence of the computed NLO cross section 
on two scales: the renormalization scale /ir and the fac- 
torization scale fip. One often makes a standard choice 
for these scales: ^ = /j,p = Pt/2. We will take this 
choice as our central value and define 

We compute the cross section near x± = X2 = 0, that is 
near the scale choice /iR = /Uf — -Pr/2. Then {xi,^} 
measures (logarithmically) the distance from this central 
value. We then fit the cross section to a quadratic poly- 
nomial P(x) in af-space, 



da{xi,X2) 




'<2cr(0,0)" 




dP T 


NLO 


dP T 


NLO 



[l + P(x)] , (3) 



where 

P(x) = X J A J + X J M JKXK , (4) 
J J,K 

with x = (xi, x%) and J, K — {1, 2}. 

We know that if we had an NNLO calculation, the de- 
pendence of the cross section on x would be canceled to 
order o? a . Thus the coefficients Aj and Mjk carry infor- 
mation about the perturbative coefficients beyond NLO. 
For this reason, we use the coefficients Aj and Mjk to 
provide an estimate of the error induced by truncating 
the perturbative expansion at one-loop order. We define 
a simple recipe for this purpose. We define an estimated 
error 3 £ sca i c as the root-mean-square average of P{x) over 
a circle with a certain radius \x\, 

4ale = ^ / d6 P{\x\ cos 9, \x\ sinf?) 2 . (5) 
27r Jo 

We need to select a value of \x\, and we make the choice 

|a?| = 2 . (6) 

In the most common method of estimating errors from 
scale variation, we would vary (2/iR/P^, 2/if/Pt) be- 
tween (1, 1) and (2, 2) and between (1, 1) and (1/2, 1/2). 
This amounts to changing x from to a vector of length 



3 We shall use £ sca i e to denote the theoretical systematic error due 
to scale dependence only, and E (no subscript) to denote the total 
theoretical systematic error. 
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\x\ = y/2 in a particular direction that corresponds to 
something close to the direction of strongest variation. 
The choice \x\ — 2 is somewhat larger than this standard 
choice. For instance, |x| = 2 in the direction x cx (1, 1) 
corresponds to 

fe^) =(2 ^^ )Sa(2 - 7 ' 2 - 7) ' (7) 

We average over the directions of x instead of taking a 
particular direction. For this reason, the value of Eq. Q 
gives results that are similar to the method that is often 
used. While varying the /i-scales along the (1,1) direc- 
tion will often work, our averaging technique provides a 
general method that seems sensible even when the one of 
the directions of slowest variation happens to align with 
the (1,1) direction. 

A straightforward calculation shows that, with the def- 
inition (|5]i, 

4.c = fi 2 +f [(TrM) 2 + 2TrM 2 ] . (8) 

We determine the coefficients Aj and Mjk by calculating 
the one jet inclusive cross section for a given value of Pt 
and rapidity. We use nine points in x-space, obtained by 
setting each {/Li R ,/xp} scale to {\Pt, \Pt, Pt} an d fit 
the results to the form given in Eq. (|3) and Eq. @. 

B. Contour plots 

We illustrate this procedure for estimating the theo- 
retical error from this source in Fig. [2] where we display 
contour plots of 1 + P(x) corresponding to the jet cross 
section at the Tevatron with Pt — 100 GeV for y = 
and for y = 2. For both values of y, we find a saddle 
point in the vicinity of {£1,22} = {0,0} which corre- 
sponds to {nn,/j,p} = {Pt/2, Pt/2}. This location of 
the saddle point is a general feature that holds through- 
out much of the kinematic range; it motivates the choice 
{/Ur,(Uf} = {Pt/2, Pt/2} as our central values. 

The estimated scale dependence error, £ sca ie> is then 
obtained by averaging the deviation of the cross section 
at a given radius in x-space. As discussed above, we 
choose a radius of |x| = 2, as indicated by the circle in 
Figure [2] The slope of the {£1,2:2} surface is steeper 
for the y = 2 case as compared with the y — case. 
Consequently, we find a larger £ sca i c for y = 2 (~ 18%) 
as compared to y = (~ 9%). 

C. Comment on the range of scale choices 

In the above analysis, we estimate the theoretical un- 
certainty by varying the fj, scales by a factor about a 
central value. This is a conventional choice, but is it 
reasonable? To examine this question, one can look at 
cases in which NNLO calculations exist. Here, we choose 
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Figure 2: Contour plot of the jet cross section in the {asi, X2} 
plane for the Tevatron (s/s = 1960 GeV) with P T = 100 GeV 
and a) central rapidity y = and b) forward rapidity y = 
2. We plot the ratio of the cross section compared to the 
central value at {xi,X2} = {0, 0}. Contour lines are drawn at 
intervals of 0.10. The (red) circle is at radius \x\ = 2. 



one typical case as an example. In Fig. [3] we show the 
NNLO cross section for Higgs production at the Large 
Hadron Collider (LHC) as a function of the P™ to pa- 
rameter as calculated by Ref. [2]. Here, the renormaliza- 
tion and factorization scales are varied by a factor of two, 
{/xr,Mf} € [A4/2, 2M h ]. 

Consider, for example, P™ to near 80 GeV. To simplify 
our argument, let us suppose that the exact QCD result 
is known and that it lies in the middle of the NNLO er- 
ror band. We then ask whether the estimated NLO error 
band was reasonable, now that we know the exact an- 



5 



pp-»H + X-»W + X-e + i/e u + X 



MRST2001 LO, MRST2004 NLO/NNLO 




P ; eto [GeV] 



Figure 3: The cross section for Higgs production at the LHC 
for LO, NLO, and NNLO calculations as taken from Ref. 
The computed cross section vetos jets (P^ > Py eto ) in the 
central region \r)\ < 2.5. 



swer. To do a real statistical analysis, we should have at 
hand many NLO calculations of separate and indepen- 
dent quantities, each with its error estimate. For each 
such quantity, a NNLO calculation that we can regard as 
nearly "exact" should be available. We would then plot 
the distribution of the differences between the NLO cen- 
tral value and the true answer in units of the NLO 1 a 
error estimate. If the error estimates are reliable, this dis- 
tribution should be a Gaussian distribution with width 
1. We cannot do that with just one datum. However, we 
can say that if the NLO estimate is reasonable then the 
central NNLO value in the one case that we have should 
be roughly 1 a away from the NLO central value. If it 
is 3 a away, then it seems likely that the NLO error was 
underestimated. If it is 0.1 cr away, then seems likely that 
the NLO error was overestimated. In the case at hand, 
the difference is about 1 cr, so we have some evidence that 
the error was correctly estimated. 



D. Scale dependence total uncertainty 

Implementing the procedure outlined above, we find 
the theoretical systematic error estimated from scale de- 
pendence, £ S caie; this is displayed in Fig. [1] The (blue) 
points are £ sca \ c computed as described above from the 
NLO cross section pQ and the (red) curve is a smooth fit 
to these points. 

We see that £ sca i e ( Pt j y), is a slowly rising function 
of Pt- For the rapidity y = at the Tevatron (y/s = 
1960 GeV), we find that £ sca le(-Pr) y) varies from 9% to 
11%. For y = 1, the uncertainty ranges from 9% to 
20%, and for y = 2 the uncertainty increases even more, 
ranging from 12% to 25% over a more limited Pt range. 



E. Scale dependence correlated uncertainty 

As described in Section [TTJ we decompose the total 
scale dependence uncertainty, £ sca i c , into a (small) num- 
ber of functions fj(Pr,y) which then combine to form 
the total uncertainty £ sca i e - 

Since the fj{Pr,y) functions represent independent 
sources of uncertainty, £ sca i c is the quadrature sum 



(9) 



We chose a set of functions Jj{PtiV) that satisfies 
Eq. ([9]). We take the fj(Pr,y) to depend on y and on 
the ratio of Pt to the quantity 4 



M(y) = y/s e- y . 
For the set of fj(Pr,y) functions we choose 



(10) 



h{P T ,y) 



h{Pr,y) 



9.62 x 10~ 2 
log(M(y)/F T ) ' 
2.89 x 10~ 2 y 2 
\og(M(y)/P T ) ' 



f 3 (P T ,y)=&A2x 10" 
U(P T , y) = 0.842 x 10- V 

MPt,ij) = 1.68 x 10- 2 log 



(11) 



/15P r \ 



\M(y)J 



f 6 (P T , y) = 0.336 x 10- V log 



( 15 P T \ 
\M(y)J 



These functions are illustrated in Fig. [5] The first two 
terms are singular as Pt — >• M(y). The first controls the 
singular behavior near y — while the second modifies 
the singular behavior for large y. The remaining terms 
constitute a polynomial in log(PV) and y 2 . Thus, we 
parameterize the y-dependence with the set of functions 
{1, y 2 }, and the Pr-dependence with the set of functions 
{1/L, 1, L} where L represents a logarithmic function of 
Pt- We believe that the parameterization in terms of 
these 2x3 = 6 functions is sufficient to reasonably de- 
scribe the theoretical uncertainties. 

Note that the coefficients of fa and f± are in the ratio 
10:1 and the coefficients of /s and fe are in the ratio 5:1. 
While we could find an excellent fit without fn and fe, 
we retain these terms to provide flexibility when one tries 
to fit the Aj coefficients to actual data. 

We can perform a similar exercise for the LHC as well; 
these results will be compiled and presented in the Sec- 
tion ELD 



4 We scale Pt by M(y) to make the argument of the logarithms 
dimensionless. This quantity provides a simple scaling, and 
roughly corresponds to scaling by the maximum Pt, P™ ax ~ 
v / 3/(2 cosh(j/)), for large y. 
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P T (GeV) P T (GeV) P T (GeV) 



Figure 4: The estimate of the uncertainty B{Pr,y) = f scale due to the scale variation as given in Eq. |8| for the Tevatron 
(y/s = 1960 GeV) with y = {0, 1, 2}. The calculation from the jet code is represented by the (blue) points, and the fit based on 
Eq. H is shown with the solid (red) curve. 




P T (GeV) P T (GeV) P T (GeV) 



Figure 5: The estimate of the uncertainty £ SC aic due to the scale variation as given in Eq. Q for the Tevatron (yfs = 1960 GeV) 
with y = {0,1,2}. The combined uncertainty £ sca i e is shown as the upper thick (red) curve, and the individual functions 
fj {Pt , y) are indicated below. 



IV. SUMMATION OF THRESHOLD LOGS 

For parton-parton scattering near the threshold for the 
production of a jet with a given Pt, there is restricted 
phase space for real gluon emission. Thus, there is an 
incomplete cancellation of infrared divergences between 
real and virtual graphs, resulting in large logarithms L 
inside the integration over parton momentum fractions. 
At n-th order in a s these logarithms enter the cross sec- 
tion in the general form a™ L 2n . The leading logarithms 
can be summed to all orders in a s . We make use of the 
numerical results from Ref. [4 , which has been imple- 
mented in the FastNLO program [S]- 

Fig. [6] displays the size of the threshold correction for 
Tevatron jet measurements at y = 0. The curve is pre- 
sented for the scale choice /i = Pt/2; we note that for this 
scale choice, the threshold correction is generally smaller 
than with other scale choices. 5 



We find the threshold corrections in this kinematic 
regime to be less than those discussed in the previous 
section (Sec. Ill I and shown in Figure [5j As the thresh- 
old corrections also arise from uncomputed higher-order 
terms, these corrections are, in a sense, already accommo- 
dated by the larger uncertainty that we estimated from 
scale variation in Eq. (111. Indeed, the functions fj for 
J = 1 and J = 2 contain singularities for Pt — > M{y) 
that are meant to incorporate the threshold singulari- 
ties. For this reason, we will not add a separate Jj(Pt, y) 
function in the expression for the total uncertainty £ to 
represent the effects of threshold logarithms. 



200 GeV, even though decreasing Pt puts us farther from the 
threshold. This rise is more pronounced for large y than we see 
for y = in Fig. [6] We suspect that this behavior is an arti- 
fact of kinematic choices in the algorithm for summing threshold 
logarithms, rather than being a real physical effect. 



5 We do not present curves for y = 1 and y = 2 because 
these curves show a rise of the correction as Pt decreases from 
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area of the jet cone — hence the factor of R 2 in Eq. ( 12 1 
At Tevatron energies, Ref. [5] finds 



A UE (1960 GeV) «3±1 GeV 



(13) 



Thus, the (Pt) shift from the underlying event correc- 
tions is given by 



(SP T ) VE « +0.7 GeV ±0.3 GeV , 
for a jet cone with R = 0.7. 

B. Hadronization correction (HC) 



(14) 



Figure 6: The ratio of the two-loop threshold resummation 
contributions for jet production compared to the total NLO 
cross section cr rcsum /(jNLO at the Tevatron (y/s = 1960 GeV) 
vs. Pt in GeV. We have set the scales to ^tp = /ir = Pt/2, 
and used y — 0. The points are computed using the imple- 
mentation of the 2-loop threshold resummation by Kidonakis 
and Owens [3]. 



The R dependence of hadronization correction is very 
different from that of the underlying event correction [6j- 
8 . The smaller the jet cone is, the more likely it is that 
hadronization will spray hadrons out of the cone. Hence, 
we will parameterize these corrections as proportional to 
1/R. Following Ref. j6], we write the hadronization cor- 
rection as 



V. UNDERLYING EVENT AND 
HADRONIZATION 

A separate source of uncertainties in jet measurements 
comes from what is colloquially known as "splash-in" and 
"splash-out" corrections. "Splash-in" corrections arise 
from the underlying event, which can deposit additional 
energy into the jet cone; we will refer to these more for- 
mally as underlying event (UE) corrections. "Splash-out" 
corrections come from the hadronization process of the 
jet which may move some of the jet energy outside the 
defined jet cone. We will refer to these as hadronization 
corrections (HC). 

In either case, the correction is modeled as adding an 
amount 8Pt to the observed transverse momentum (or 
transverse energy) of the jet. We denote the average over 
many events of 5Pt by (SPt)- A complete analysis of the 
UE and HC contributions was performed by Cacciari, 
Dasgupta, Magnea, Salam in Refs. [SHE]. We find this to 
be an entirely suitable method for our estimate of (SPr), 
and we adapt their results in the following. 



A. Underlying event (UE) 

We can parameterize the effect of the underlying event 
corrections on the apparent Pt of the jet as 



(SPt) ue = Mje^R 2 



(12) 



where R is the cone radius of the jet and Afje is the 
average transverse energy per unit rapidity in the under- 
lying event. Because we model the "splash-in" energy as 
random and uncorrelated with how the jet develops, the 
contribution from the underlying event will scale as the 



(SB 



T/HC 



(15) 



where -4(/i/) parameterizes the soft gluon radiation. 
Ref. takes /x 7 = 2 GeV, and finds A(2 GeV) « 
0.2 GeV. In Eq. (15 1, Cj is a color factor that depends 
on whether the jet is initiated by a quark, for which 
Ci = Cp =4/3, or by a gluon, for which C L = Ca = 3. 
We thus need an estimate of the fraction of jets that are 
gluon jets. Using calculations from the literature [5], we 
estimate that, for the Tevatron in the low Pt region, the 
fractions of quark and gluon jets are approximately 



fa 



Using these fractions, we can form a weighted average of 
the quark and gluon terms to obtain 



(SP T )uc = f g (SP^uc + f g (SP^nc 

f 2Cf AC \ t 2Ca 
= -Jq—j^A(pi) - fg 

-1 GeV ±0.5 GeV 



R 



A(m) (16) 



Here, we have used a typical cone radius of R = 0.7 and 
taken a conservative choice for the uncertainty of 50% of 
the correction. 



C. (SP T ) from the UE and HC 



Combining the underlying event of Eq. ( 14 1 and the 
hadronization corrections of Eq. (16), the net Pt shift is 



(5Pt 



-0.3 GeV ±0.6 GeV 



(17) 



where we have added the separate uncertainties in 
quadrature. 
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TOT 




Figure 7: We display the expected Pt shift, (SPt), in GeV 
vs. jet cone radius R for the UE, HC, and combined results 
(TOT) at the Tevatron. The calculation of the HC uses a 
combination of quark-initiated ( f q — 2/3) and gluon- initiated 
(f g = 1/3) jets. The upper solid (blue) line represents the UE 
correction, and the lower solid (green) line represents the HC 
terms. The combination of these corrections (TOT) is repre- 
sented by the central (red) band including the uncertainties. 
The vertical line corresponds to R=0.7. 
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Figure 8: Jet cross section cPa/dPr/dy vs. Pt in GeV with 
y — at the Tevatron in units of nb/GeV. The line is a power 
law fit with n — 7; this describes the slope of the jet data in 
the range P T « [50, 300] GeV. 



The individual underlying event and hadronization re- 
sults for (SPt) are displayed in Fig. ff\ fo r the Tevatron 
using the parameterizations of Eq. ( |14[ ) and Eq. (16 1. 
The combined result for (SPt), including the uncer- 
tainty band, is also displayed. The underlying event and 



hadronization corrections have opposite sign, and we note 
that for a jet cone radius of R = 0.7, the two corrections 
nearly cancel each other. 



D. From (SP T ) to 5a 



dcr(P T ) 
dP T 



const. 
pn 



(18) 



in the specific Pt range of interest. For jets at the Teva- 
tron in the intermediate Pt range of ~ [50, 300] GeV, we 
find n ~ 7 as illustrated by Fig. [8] 

The effect of the underlying event and hadronization 
corrections is to shift the jet Pt from its value P^ ert at 
the NLO parton level to a new value 



P T 



ppcrt 
T 



(SPt 



where (SPt) is the average change in the transverse jet 
transverse momentum due to underlying event additions 



and hadronization subtractions from Eq. ( 17 1 



If we write the true differential cross section as a func- 
tion /, 

dP T ~ n ^ T) ' 

then / is related to the perturbatively calculated function 
/pert by 

f(P T ) « / P crt(P| crt ) = /p-t (Pt - (SPt)) ■ 

We can perform a Taylor expansion about Pt for small 
SPt, 



/(Pt) « f P crt(P T - (SPt)) 



/pert(^r) - (SPt 

/pert (P T )<\ 



dP T 
(SPt) ] 
Pt J ' 



Here we have used the power law of Eq. ( 18 1 to replace 
f(P T ) by -nf(P T )/Pr- Thus, to first order we find 6 



da 
d~Pr 



da 



pert 



dP 7 



1 



(SPt 
Pt 



(19) 



so that the fractional correction is n (SPt) /Pt- Using 
n fa 7 and the estimate from Eq. (17) of (SPt), we find 
that the fractional correction to the cross section is ap- 
proximately 



7 x 



-0.3 GeV ±0.6 GeV 
Pt 



2 GeV 4 GeV 



Pi 



Pi 



Thus we estimate the fractional uncertainty from the un- 
derlying event and hadronization to be 4GeV/Pr. 

We account for this source of uncertainty by adding a 
new function fj(Pr,y) with J = 7, 



h(P T ,v) 



4 GeV 
Pt 



(20) 



for Tevatron jets in the Pt range of ~ [50, 300] GeV. 



The differential jet cross section can be approximated 
by a power law of the form 



6 Cf., Eq. (5.9) of Dasgupta et al. in Ref. [6] 
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Uncertainty fj 


Source 


{.h, h, /3, fi, fs, h} 


perturbative 


h 


non-perturbative 



Table I: A compilation of the source of uncertainties (Jj) that 
comprise the total jet cross section uncertainty £. The per- 
turbative uncertainties arise from the higher, uncalculated, 
orders of perturbation theory and are estimated using the 
{/ZF,/iR.} scale variation of the calculated cross section. The 
non-perturbative uncertainties are an estimate of the under- 
lying event and hadronization corrections. 



VI. SUMMARY FOR THE TEVATRON 

We have described the correlated theoretical system- 
atic uncertainty using a total of seven functions, as sum- 
marized in Table [I] The net error at any one value of 
{Pr,y} is obtained by adding these seven functions in 
quadrature 



£(P T ,y) = JJ2fj(P T , y y 



(21) 



We now summarize the complete set of contributions 
to the uncertainty of the differential jet cross section as 
a function of {Pt, y} for the Tevatron: 



h{PT,y) 



h{P T ,v) 



9.62 x 1CT 2 
log(M(y)/P T ) ' 
2.89 x lCT 2 y 2 
\og{M{y)/P T ) ' 



f 3 (P T ,y) = SA2x 10 

U(P T ,y) 



0.842 x I0~ 2 y 2 



h(P T ,y) 
h{P T ,y) 
fr(P T ,y) 



L68xl0 ~ 2 log (H) 



(22) 



0.336 x 10" 
4GeV 



y 2 log 



( 15Pt \ 
\M(y)J 



Pt 



We display these results in Figure [9 For Pt > 100 GeV, 
the perturbative uncertainties arc dominant, and slowly 
rise with increasing Pt; this results holds across the full 
y-range, but the rise with Pt is more pronounced at large 
y. For P T < 100 GeV, the uncertainty from the UE 
and HC terms become increasingly important as Pt de- 
creases. 



VII. THEORY ERRORS AT THE LHC 

Having demonstrated the method for determining the 
theoretical systematic uncertainty at the Tevatron, we 
perform a parallel analysis for the Large Hadron Collider 
(LHC). 



A. Perturbative uncertainty 

We again estimate the error from not having calculated 
beyond NLO by using the dependence of the NLO cross 
section on the scales {/XR,/ip}, just as in the Tevatron 
case, and this yields the functions {f 1} f 2 , f 3 , f 4 , fs, fe} 
summarized in Eq. d27| at the end of this section. 



B. Underlying event and hadronization 

We proceed as in Sec. [V] for the Tevatron, accounting 
for the changed circumstances at the LHC. We first need 
to estimate the error in the determination of the contri- 
bution to the average jet transverse momentum, (SPt), 
arising from the underlying event and from hadroniza- 
tion. 

The underlying event contribution to (SPt) is deter- 
mined by the parameter Aue in Eq. ( fl2| ). Consistently 
with Refs. [SHE], for the LHC we take~A UE (14 TeV) 
10 ± 4 GeV, and obtain 



-2.5 GeV ± 1 GeV 



(23) 



For the contribution to (SPt) from hadronization, we 
use Eq. (16 1 with A(^i) ~ 0.2 GeV as before. For the 
fractions f q and f„ of quark and gluon jets in the rela- 
tively low Pt region where the hadronization corrections 
are significant, we use 



Using these fractions, we can form a weighted average of 
the quark and gluon terms and estimate the hadroniza- 
tion contribution to (SPt) to be 



(SP T ) HC = -1.4 GeV ±0.7 GeV 



(24) 



Combining the underlying event and hadronization 
contributions, we estimate 



(SPt 



4 GeV ± 1.2 GeV 



(25) 



where we have added the separate uncertainties in 
quadrature. 

The results for the underlying event and hadronization 
contribution to (SPt) are displayed in Fig. 10 for the LHC 



using the parameterizations of Eq. ( 24 1 and Eq. ( 23 ) but 
with a variable cone size R. 

The correction to (SPt) determines the correction to 
the cross section via Eq. (19 1. For this, we need the 



power n that describes the approximate power law fall- 
off of the cross section. As illustrated in Fig. [TT] a power 
law with n ss 6 describes the data over the range Pt ~ 
[100, 1000] GeV. Using n « 6 and the estimate from 
Eq. (25) of (SPt), we find that the fractional correction 



to the cross section is approximately 
1 GeV ± 1.2 GeV 



6 x 



Pi 



6 GeV 7 GeV 



10 




P T (GeV) P T (GeV) P T (GeV) 

Figure 9: A compilation of the uncertainties for jet production at the the Tevatron (^fs = 1960 GeV) for y = {0,1,2}. The 
numeric label corresponds to the error components summarized in Eq. ( |27[ ). The upper thick (red) line is the quadrature sum 
of the individual errors . 



TOT 




10- 



Figure 10: We display the expected Pt shift, (SPt), in GeV 
vs. jet cone radius R for the UE, HC, and combined re- 
sults (TOT) at the LHC. The calculation of the HC uses a 
combination of quark-initiated (f q — 1 /3) and gluon-initiated 
(f g = 2/3) jets. The upper solid (blue) line represents the UE 
correction, and the lower solid (green) line represents the HC 
terms. The combination of these corrections (TOT) is repre- 
sented by the central (red) band including the uncertainties. 
The vertical line corresponds to R=0.7. 



Thus we estimate the fractional uncertainty from the un- 
derlying event and hadronization to be IGeY /Pt- We 
include this in the estimate of systematic theoretical er- 
rors by including a function Jt{Pt) given by 



h{P T ) 



7 GeV 
Pt 



(26) 



for LHC jets in the range P T w [100, 1000] GeV. 



C. Summary: LHC 

We now summarize the complete set of contributions 
to the uncertainty of the differential jet cross section as 



C 
O 



y 10 



u 

X/i 



- 2 



GO 

2 io- 4 
u 

5 IO" 6 



10" 



LHC 
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Figure 11: Jet cross section d 2 a / cIPt / dy vs. Pt in GeV with 
y = at the LHC {s/s = 14 TeV) in units of nb/GeV. The 
line is a power law fit with n = 6; this describes the slope of 
the jet data in the range P T « [100, 1000] GeV. 



a function of {Pt, y} for the LHC: 
4.56 x 10~ 2 



h(PT,y) 



log(M(y)/P T ) 
1.24 x 10~ 2 y 2 
log(M(y)/P T ) 



/ 3 (Pt,2/) = 5.36x 10 
U{Pt,v) 



0.536 x 10~ V 



f 5 (P T ,y) = LOT \ 10 
h{P T ,y) 
h{P T ,y) 



-2 1 f 15P A 

0.214 x 10- 2 y 2 logY^^ 



(27) 



7 GeV 



Pi 



We display these results in Figure 12 In the central ra- 
pidity (y ~ 0) region for Pt > 500 GeV the perturbative 
uncertainties are dominant and slowly rise with increas- 
ing Pt, while for Pt < 500 GeV the nonperturbative 
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Figure 12: A compilation of the uncertainties for jet production at the LHC (y/s = 14, 000 GeV) for y — {0, 1, 2}. The numeric 
label corresponds to the error components summarized in Eq. (27 1. The upper thick (red) line is the quadrature sum of the 
individual errors. 



uncertainties become increasingly important. For y = 2, 
the transition Pt is closer to 300 GeV than 500 GeV. 

VIII. CONCLUSIONS 

As the LHC prepares to take data, it is important that 
we be able to determine whether a physics signal is consis- 
tent with the standard model. For example, if we observe 
a signal that is inconsistent with the standard model pre- 
diction, but the augment for this inconsistency includes 
only experimental errors, we cannot claim this is "new 
physics" until we demonstrate it is also inconsistent in- 
cluding both experimental and theoretical errors. This 
paper provides a framework to quantitatively make such 
a determination in the case of jet physics. Similarly, this 
paper provides a framework to quantitatively fit parton 
distribution functions to Tevatron and LHC jet data, in- 



cluding estimated errors from the theory. 

The framework that we provide involves functions 
fj(Pr, y) that represent independent contributions to the 
theory error. We note that other authors might estimate 
the errors differently and thus produce different functions 
fj(Pr,y)- We hope that this will happen and that the 
merit of different choices will be debated. 
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